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Simple analytical relations for reflection and transmission matrices of plane-parallel layers of random media 
with discrete particles are presented. They can be used for rapid estimation of intensity and polarization 


characteristics of reflected and transmitted light beams under arbitrary illumination of a layer. 
of the analytical formulas obtained increases with the optical thickness 7 of a layer. 
Another limitation is due to the probability of photon absorption £, 


plicable only at large values of 7(7 > 5). 


The accuracy 
Thus equations are ap- 


which should be rather low (8 < 0.05-0.1, depending on the optical characteristics in question). © 2001 Op- 


tical Society of America 


OCIS codes: 010.1290, 010.1310, 120.0280, 260.2130, 290.4210, 290.7050. 


1. INTRODUCTION 


Reflection and transmission matrices of various disperse 
media are usually found by solving the integro- 
differential vector radiative transfer equation.! The ra- 
diative transfer equation can be applied to a broad range 
of physical problems, including neutron transport, atmo- 
spheric, oceanic, and tissue optics. It is solved numeri- 
cally in most of cases. However, analytical solutions of 
this equation have also been derived for many particular 
situations, e.g., for optically thin and optically thick 
layers.! Formulas for thin layers are especially simple 
and have already been used in many applications.? The 
use of asymptotic analytical equations for optically thick 
layers is not so common, because asymptotic formulas are 
expressed in terms of solutions of integral equations.! 
One can argue that there is no point in using asymptotic 
solutions if they require complex numerical codes for cal- 
culation of the auxiliary functions and parameters in- 
volved. Perhaps it is better to start from the radiative 
equation from the very beginning in this case. 

However, simple relations between global optical char- 
acteristics (e.g., the reflection function and the degree of 
polarization) of different disperse media and their local 
optical characteristics (e.g., the absorption coefficient and 
the phase matrix) are especially important for the solu- 
tion of inverse problems. Thus it is the aim of this paper 
to derive such equations for weakly absorbing plane- 
parallel optically thick light-scattering layers. The as- 
sumption of low level of light absorption in a medium is 
used to simplify the general asymptotic solutions, which 
have already been known for a long time.** Note that 
weakly absorbing optically thick layers occur frequently 
in nature. Snow and ice fields, fog and clouds, blood, pa- 
per, paints and many other types of natural media belong 
to this important subclass of disperse media. 
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2. ASYMPTOTIC THEORY 


Let us start from asymptotic formulas of the vector radia- 
tive transfer theory.! Corresponding equations can be 
written in terms of reflection R and transmission 7 ma- 
trices. These matrices relate Stokes vectors of transmit- 
ted S, and reflected S, light beams to the Stokes vector of 
the incident light beam F: 


A 


Si(u, $) = T(p, Mo, oF, 


S.(u, 6) = Ê(u, po, oF 


Here u = cos 0; mo = cos 6; 0o and 0 are incidence and 
observation angles, respectively; ¢ is the relative azimuth 
of incident and reflected or transmitted light beams; and 
we assume that light comes from just one direction speci- 
fied by the cosine of the incidence angle my and the azi- 
muth ġo. Note that the vector F is normalized in such a 
way that the first element of 7F is the net flux of a light 
beam per unit area of a disperse medium. 

Thus we will use the following asymptotic 
expressions!" for the reflection matrix R( LL, Ho, $) and 
the transmission matrix Î( Lt, ko) of a homogeneous 
plane-parallel turbid layer of large optical thickness r 
= CoxtL (Cex, is the extinction coefficient of the disperse 
slab in question, and L is the geometrical thickness of the 
slab): 


R(n, no, $) = Rul, no, $) — fT (Cu, uo). (1) 


(u, uo) = tK(u)K? (u0), (2) 

where 
m exp(—k7) 3 
t= ae pe (3) 
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It should be pointed out that matrices Ê H, Ho, $) and 
T(u, uo) can be used to calculate the Stokes vectors of re- 
flected and transmitted light beams under arbitrary illu- 
mination of a scattering layer.! Equations (1) and (2) fol- 
low from the exact vector radiative transfer equation as 
7— ©, Thus they are exact formulas in the asymptotic 
limit. The details of the derivation of these important 
equations were presented for the first time by Domke.* 
Note that the vector K7( uo) is transpose to the vector 
K(uo). Constants f and m and the vector K() can be 


found from following relations*”: 
f = s exp(—k7), (4) 
1 
m = 2| duu[P(u)P(u) — P™(-p)P(—-z)], (5) 
0 
1 A 
Ku) = m~"|P(u) - 2 I. dééR.(u, €)P(—€)|, (6) 


where P? means the transpose vector to the vector P, 
1 
s= 2| duuK”(u)P( — n), (7) 
0 


and 


A 


1 27 = 
R.(m, £) = =|, dGR..(u, € $) (8) 


is the azimuthally averaged reflection matrix of the semi- 
infinite medium with the same local optical characteris- 
tics as those for the finite layer under study. Note that 
the escape vector K(,.) also occurs in the so-called Milne 
problem. It describes the angular dependence of the in- 
tensity and polarization characteristics of light emerging 
from semi-infinite turbid layers with sources of radiation 
located deep inside the medium. The constant k is the 
smallest discrete eigenvalue that is real and nondegener- 
ate and obeys the following equation: 


wo f1 7 
(1 — kp) P(p) = 2 f ateu, EPE), (9) 


which is called the vector characteristic equation of the 
radiative transfer theory. Here y is the cosine of the ob- 
servation angle, wọ is the single scattering albedo, which 
is equal to the ratio of the scattering and extinction coef- 
ficients, Z(, £) is the azimuth-averaged phase matrix, 
and P(€) is the corresponding eigenvector. 

Equations (1) and (2) are very general and can be ap- 

plied to disperse media with arbitrary local optical char- 
acteristics, assuming that the optical thickness of the 
light-scattering medium in question is large. However, 
their application to the solution of practical problems is 
not easy, because it is necessary to solve integral equation 
(9) and perform integrations (5)—-(8). The reflection ma- 
trix R(n, é, &) obeys the integral equation as well.’ 
F The problem of finding matrices R(u, uo, 6) and 
T(u, mo) in Eqs. (1) and (2) can be simplified for a broad 
class of so-called weakly absorbing media, which are char- 
acterized by a small value of the probability of photon ab- 
sorption 8 = 1 — wọ. One can obtain,! at small values 
of b; 
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k = [8(1 - g)(1 - w)]”, (10) 
8k i 
m = 3a -a p (11) 
4ka 
$22 d= (12) 
3(1 — g) 

; a 4k _ r 
R(n, Mo, $) = Rolu, Mo, $) — z~~ Kolu)K (uo), 
—— (13) 


where K)() and Ên, Ho, $) are the escape vector and 
the reflection matrix, respectively, of the nonabsorbing 
medium with the same phase matrix as the finite absorb- 
ing layer under study, 


1 fr 

g= zl p(0)sin 0 cos 6d0 (14) 
0 

is the asymmetry parameter of the phase function p( 9), 

and 


1 
a= Ji duu? S(u), (15) 
0 


where S(u) = K3(w)e,, e; is the unit column vector. It 
should be pointed out that integral (15) is close to 1 for all 
types of phase matrices.*? Thus we will neglect small 
differences A = 1 — æ in the discussion that follows. 
Note that the first moment of the function S(u) does not 
depend on the phase matrix at all, because of the normal- 
ization condition?: 


1 
af duuS(u) = 1. (16) 
0 


Thus Eqs. (10)-(18) allow us to reduce the complexity 
of Eqs. (1) and (2) Only function K(u) and 
RX H, to, $) for a nonabsorbing medium should be found 
numerically in this simpler case. They depend on the 
phase matrix of the random medium in question. How- 
ever, they do not depend on the single-scattering albedo 
wo and the optical thickness 7. Another interesting fea- 
ture of these functions is their weak dependence on the 
microstructure parameters (size, shape, and chemical 
composition of particles) of the medium under study in 
the broad range of angular parameters u, wo, and ¢. 
This is due to the randomization of both the polarization 
states and the propagation directions of photons in highly 
scattering semi-infinite layers irrespective of the local op- 
tical properties of the medium in which they propagate. 


3. MODIFIED ASYMPTOTIC EQUATIONS 


A. Main Formulas 

Numerical calculations show that expansions (11)-(13) 
can be applied only for extremely small values of 6 = 1 
— wo (B~ 10%). This greatly limits the power of the 
asymptotic theory for real-life problems. 

Various methods to overcome this problem have been 
used. The most straightforward approach calls for the 
derivation of higher-order terms in expansions (10)—(13). 
However, here we will use another method, which was 
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initially applied in scalar radiative transfer theory. Ac- 
cording to this alternative approach, expansions (11) and 
(12) are replaced by exponential functions, 


m = 1 — exp(-2y), (17) 
s = exp(—y), (18) 
where 
4k 
y= 30 —)’ (19) 


and we assume that a= 1. Equations (17) and (18) 
should be regarded as empirical relations; they cannot be 
derived from first principles. Formulas (17) and (18) 
transform to exact asymptotic results (11) and (12) as 
k — 0, provided that a = 1. These equations do allow us 
to calculate radiative and polarization characteristics for 
many types of absorbing light-scattering media with val- 
ues of 8 < 0.1, as follows from the results presented be- 
low. 

Thus Eqs. (1) and (2) can be transformed into the fol- 
lowing forms: 


R(u, no, $) = Ês(u, uo, $) 
— Î(u, wo)exp(-y — kr), (20) 


T (1, mo) = tKo(4)KG( 0), (21) 


where we have accounted for Eqs. (17) and (18) and the 
approximate equality 


mK()K" (419) ~ [1 — exp(—2y)]Ko( 4) Kj (uo). 


(22) 
The value of 
1 — exp(—2y) 
t= (23) 
exp(kr) — exp(—2y — k7) 
is the global transmittance, defined as! 
1 1 
t= af nan | HoduoTi(u, mo). (24) 
0 0 


Equation (24) follows from Eq. (21), accounting for the in- 
tegral relation,! 


1 
af uduKo(u)du = 1, (25) 
0 


where Ko;() is the first component of the escape vector 
Koy). 

Note that the global transmittance (24) could be trans- 
formed into the following simple form®: 


shy 
t= err (26) 
where 
x = kr. (27) 


Equations (20) and (21) are much simpler than initial 
asymptotic formulas (1) and (2). They allow us to calcu- 
late the reflection and transmission matrices of optically 
thick weakly absorbing disperse media by simple means, 
if the solution of the problem for a semi-infinite medium 
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with the same phase matrix as for the finite layer under 
consideration is available. This reduction of the general 
problem to the case of a semi-infinite medium is impor- 
tant for radiative transfer theory. 

It follows from Eqs. (20) and (21), neglecting the polar- 
ization characteristics of a light beam under study, that 


R(n, Mo, $) = R(n, Mo, $) — T(u, wo)exp(-y — k7), 
(28) 


T(u, Mo) = tKo(u)Kolno), (29) 


where T = T11, R = Ry, R? = R2,,, Ko = Ko are the 
first elements of the corresponding matrices. Equations 
(28) and (29) represent well-known results of scalar radia- 
tive transfer theory.*” This confirms our calculations. 

Note that formulas (20) and (21) can be simplified for 
nonabsorbing media (k = 0): 


R(n, Ho, 6) = Êu, Ho, $) — T(w, Mo), (80) 


T(u, oo) = tKo()K5 (uo), (31) 
where [see Eq. (23) as k — 0] 
t= 141+ 271 - g)]. (32) 


Let us apply Eqs. (30) and (31) to a particular problem, 
namely, to the derivation of a relation between the spheri- 
cal albedo r = 1 — ¢ and the degree of polatization of 
reflected light P(uo) at the nadir viewing geometry 
(u = 1). We will assume that a scattering layer is illu- 
minated on the top by a wide unidirectional unpolarized 
light beam. The value of P(yuo) is given simply by 
—Rə:( uo, 1)/Ry4(49,1) in this case. Thus it follows 
from Eqs. (30) and (31) that 


P..(Mo) 


P(uo) = Cada (33) 
where 
Ko1(1)Koi( uo) Ruoi(Ho, 1) 
oo Roi1(Ho, 1) aa Ronluo, iy" 
(34) 


and we have accounted for the equality” Ko9(1) = 0. Our 
calculations show that the value of u( mọ) is close to 1 for 
most of the observation angles, which implies inverse pro- 
portionality [see Eq. (33)] between the brightness of a tur- 
bid medium and the degree of polarization of reflected 
light. 

This inverse proportionality between the spherical al- 
bedo r and the degree of polarization P()) was discov- 
ered experimentally by Umow® almost 100 years ago. 
Equation (33) can be considered a manifestation of this 
important law, which has important applications in plan- 
etary spectroscopy.” 


B. Numerical Results 

Let us apply the equations derived to the specific case of 
water clouds illuminated by solar light. Numerical cal- 
culations with Eq. (34) show that the function u( uo) var- 
ies from 1 to 1.36 with the incidence angle in this case. 
However, this variation is of secondary importance in cal- 
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culations of the degree of polarization because of small 
values of the difference 1 — r for weakly absorbing opti- 
cally thick media [see Eq. (33)]. Thus one can neglect the 
angular dependence u(y.) and use the average value 
u( uo) = 1.18 in the case under study. It follows that Eq. 
(33) can be written in the following simplified form: 


P..( Lo) 
P(o) = TEET (35) 


where global transmittance ¢ is given by Eq. (32). Let us 
study the accuracy of this simple equation. The accuracy 
of Eq. (35) can be seen from Fig. 1, where data obtained 
from this simple formula and results of the numerical so- 
lution of the vector radiative transfer equation for light- 
scattering media with water droplets by the doubling 
method! are presented. The value of P..(u9) was ob- 
tained from the exact radiative transfer calculations for a 
semi-infinite layer. The particle size distribution (a) in 
a cloud layer was given by the following formula: 


a 
-u — 
ao 


w(a) = Aa“ exp 


? 


where A is the normalization constant (Jof(a)da = 1), a 
is the radius of particles, ag = 4 um, and u = 6. This 
distribution is often used in cloud optics. It is called the 
gamma size distribution. The optical thickness of a layer 
was equal to 8, 15, and 30. Calculations of the phase ma- 
trix were performed with Mie theory at a wavelength of 
443 nm, where light absorption by water droplets can be 
neglected. The obtained value of the asymmetry param- 
eter was equal to 0.8541 in the case under investigation. 

It follows from Fig. 1 that Eq. (85) describes the degree 
of polarization of reflected light with high accuracy at 7 
= 8. Asa matter of fact, this equation can be applied to 
media with optical thickness as low as 5 (see Fig. 2). 

Equations (33) and (34) can be easily generalized to ac- 
count for the absorption of light in a medium: It follows 
from Eq. (20) for the value of P = —Ro1(puo,1)/ 
Rl uo, 1) that 


R.221(Ho, 1) 
Rs11(u0, 1) — Ko(u)Koi(uo)t exp(—x — y) 


P(o) 


or 


PŽ( uo) 


P _ 
oa 1 — w*(o)t exp(—x — y)’ 


(36) 


where the global transmittance ¢ is given by Eq. (26) and 


(440) Ko(1)Ko1( uo) (37) 
u ee 
Ao R%11(40, 1) 


Values of P%ž(uo) and Rž(uo, 1) represent the degree of 
polarization and the reflection function of a semi-infinite 
weakly absorbing medium at the nadir viewing geometry. 
Note that Eq. (36) can be written in the following form: 


P(uo) = €( Mo, T)PE(Mo); (38) 


where the value of 
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1 
,T) = 39 
eC 1 — u*(mo)t exp(—x — y) = 


can be interpreted as the enhancement factor, which is 
due to the finite depth of a turbid layer. Equation (38) 
can be used to estimate the change in the degree of polar- 
ization with the optical thickness of the medium in ques- 
tion. 

It is evident that the value of c — 1 as T— œ, as it 
should be. Also it follows from Eq. (38) that zeros of po- 
larization curves for semi-infinite and optically thick fi- 
nite layers coincide, which is supported by numerical cal- 
culations with the radiative transfer code (see Fig. 1). 
This is due to the fact that multiple light scattering fails 


20 ere eee : , gerien 


15 F p z 


degree of polarization, % 


-5 f 1 1 1 j 
0 10 20 30 40 50 60 70 80 90 


incidence angle, degrees 


Fig. 1. Dependence of the degree of light polarization reflected 
from a water cloud with gamma particle size distribution at a 
wavelength of 443 nm on the incidence angle at nadir observa- 
tion at 7 = 8, 15, and 30 according to numerical radiative trans- 
fer calculations (broken curves) and Eqs. (35) and (32) (symbols). 


30 EIEEE T r =e 


approximate 
© exact 


20 + oO 4 


degree of polarization, % 
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0 1 1 : , f 
0.00 0.05 0.10 0.15 0.20 0.25 
inverse optical thickness 


Fig. 2. Dependence of the degree of light polarization reflected 
from a water cloud with gamma particle size distribution at a 
wavelength of 443 nm on the inverse optical thickness at nadir 
observation and incidence angle equal to 37° according to ap- 
proximate Eq. (35) (solid curve) and numerical radiative transfer 
calculations (symbols). 
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Fig. 3. Dependence of the degree of light polarization reflected 
from a water cloud on the probability of photon absorption at na- 
dir observation and incidence angle equal to 30° and 37° accord- 
ing to numerical radiative transfer calculations (symbols) and 
Eq. (86) (lines) at optical thickness 7. The phase matrix coin- 
cides with the phase matrix used for calculations presented in 
Figs. 1 and 2. 


to polarize incident unpolarized light. It only diminishes 
the polarization of singly scattered light. Thus the 
angles where the polarization is equal to zero for singly 
scattered light are preserved for semi-infinite layers as 
well. 

The dependence of the degree of polarization at fixed 
incidence angles as a function of the probability of photon 
absorption, calculated with Eq. (36) and the vector radia- 
tive transfer code, is presented in Fig. 3 at 4 
= 30°,37°,0° for the same phase matrix as in Figs. 1 and 
2. One can see that the accuracy of Eq. (36) is better 
than 2% at 8 < 0.1 at these geometries. 


4. CONCLUSIONS 


In conclusion, Eqs. (20) and (21) can be used for studies of 
the reflection and transmission matrices of optically thick 
disperse media, provided that the reflection matrix of a 
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semi-infinite medium with the same microstructure is 
known. The phase matrix of the scattering medium un- 
der investigation can be arbitrary. However, the absorp- 
tion of light by particles should be rather low, which is 
usually the case in the visible region of the electromag- 
netic spectrum. 

The formulas obtained open new possibilities for the 
parameterization of global optical characteristics of dis- 
perse media. They also could be useful in the solution of 
inverse problems, particularly those that involve spectro- 
polarimetric measurements. 
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